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I. INTRODUCTION 

One of the fundamental models in quantum optics is 
the interaction of a spin- 1/2 system with a bosonic mode 
Q . The most prominent example is cavity quantum elec- 
trodynamics (CQED), where a two level atom interacts 
with a single mode of the radiation field in a high-Q cav- 
ity. CQED has been the topic of a series of seminal exper- 
iments both in the microwave and optical regime, demon- 
strating quantum control on the level of single atoms and 
photons in an open quantum system 0, 0i El- • 

In the present paper we will consider a system with 
the same basic ingredients, however in the context of cold 
atoms and quantum degenerate gases. The key feature of 
these systems is there controllability and weak decoher- 
ence. In particular we employ two aspects of control, the 
confinement of atoms in optical lattices [E 0> H> El an d 
(magnetic or optical) Fcshbach resonances as a way to 
manipulate atomic interactions |ToL ITU Il2l Il3| . Accord- 
ing to the setup described in Fig. OJa) we will study the 
dynamics of an atomic quantum gas in ID (with a single 
internal atomic state), representing bosonic or fermionic 
"modes" , controlled by an atomic spin-1/2 impurity. The 
quantum gas is confined by tight trapping potentials (e.g. 
an optical or magnetic trap), so that only the motional 
degrees along the z-axis in Fig.^a) are relevant. In the z- 
direction the motion is confined to the left by a trapping 
potential (e.g. a blue sheet of light), while the atomic 
impurity restricts the motion of the gas to the right due 
to collisional interactions of the quantum gas with the 
impurity. The atom representing the impurity can, for 
example, be a different atomic species in a tight trapping 
potential, a configuration discussed in Refs. [Til llq as 
an atomic quantum dot (OD system). Thus the impurity 
atom plays the role of "single atom mirror" confining the 
quantum gas in an "atomic cavity" . 

In our model system the impurity atom is an internal 
two level system, which we write as an effective spin- 
1/2. In the following we will also interpret this two- 
level system as a qubit with two logical states |0) = | [) 
and |1) = I |). Cold atom collision physics allows for 
a situation where the collisional properties (scattering 




FIG. 1: (a) A spin 1/2 impurity Q used as a switch: in one 
spin state it is transparent to the probe atoms A, but in the 
other it acts as a single atom mirror, (b) Implementation 
of the SAT: The impurity atom Q and the probe atoms A 
are trapped separately in state-dependent ID optical lattices. 
The probe atoms A initially are in a Mott insulating state. 



length) of the impurity atom and atoms in the quantum 
gas are spin-dependent. As illustrated in Fig. |TJa), we 
assume that in one spin state, say | |), the single impurity 
atom is completely transparent for the quantum gas, i.e. 
the gas will leak out through the "mirror" . In contrast, in 
the other spin state the mirror atom is "highly reflective" 
confining the gas. For an impurity atom (qubit) initially 
prepared in a spin superposition 

|Vo(t = 0)) = a t | D+aJ |) 

the combined system at a time t will be in a macroscopic 
superposition state 

|*(t)> = a T |T>|<M*)> + "lU>l<M*)>' (i) 

with \<f> a (t)) many body wave functions of the gas atoms. 
Thus \^{t)) represents a Schrodinger cat state of two en- 
tangled quantum phases of gas atoms, the first one corre- 
sponding to gas confined by the mirror (Fig. ^a) upper 
figure) and the second one to the expanding gas (Fig.QJa) 
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FIG. 2: The single atom mirror as a macroscopic quantum 
gate. The qubit Q entangles two distinguishable macroscopic 
phase of the probe atoms A and provides for a quantum non- 
demolition interaction. 



lower figure) . The entanglement of the spin with a macro- 
scopic number of atoms can be interpreted as a macro- 
scopic quantum gate, as explained in Fig. EJ) , implement- 
ing a quantum nondemolition interaction (QND) |16| . 
In this sense the setup represents a "amplifier" of the 
state of the qubit. This situation is reminiscent of a 
Single Electron Transistor (SET) in mesoscopic physics 
|17|. and has stimulated the name Single Atom Transis- 
tor (SAT) for the setup Fig. Ufa) in Ref. 0, with the 
essential difference that the dynamics underlying Q is 
completely coherent. We finally remark that this setup 
also allows for a single shot QND measurement of the im- 
purity atom (qubit) by observing in a single experiment 
the distinct properties of the |</>j.(i)) or (^-(t)} quantum 
phases. 

As a variant of the configuration of Fig. da) we will 
consider below in particular the case where the quan- 
tum gas is loaded in an optical lattice, as illustrated in 
Fig. Gib) . In this case the gas could be loaded initially, 
for example, in a Mott insulating state, i.e. where large 
repulsion of the gas atom leads to a filling of the lattice 
sites with exactly one atom per lattice site |7|, 13, [l9j . The 
cat state will thus correspond to a superposition of 
the Mott phase and the melted Mott phase, i.e. a (quasi-) 
condensate of atoms obtained by expansion of the atomic 
gas: 



|*(t))=a i U)|BEC)+a T |T)|Mott). 



(2) 



In this case the distinguishing features of the two en- 
tangled quantum phases are the observation / non- 
observation of interference fringes as signatures of the 
Mott and BEC phase, when the atomic gas in released 
in a single experiment. 

Transport through an impurity is a well studied prob- 
lem in mesoscopic condensed matter physics [2(J, |2j, |23, 
l23l ] , which typically focuses on conductance properties of 
a system attached to leads. In contrast, in the context 
of cold gases we have a full time-dependent coherent dy- 
namics in an otherwise closed system. 

A short summary of the present work including results 
from numerical studies was presented in Ref. [l8| . In 
this paper we will present details of our analytical calcu- 
lations, while we refer to Ref. p^j on a complementary 



numerical treatment of these problems using time depen- 
dent DMRG techniques. 

The paper is organized as follows: In Sec.^we intro- 
duce the model used for describing the implementation 
of the Single Atom Mirror using cold atoms in optical 
lattice. In Sec. IIIII we consider the detailed scattered 
processes involved in the transport of a single particle 
through the mirror. We solve exactly the scattering prob- 
lem in the lattice by integrating the Lippmann-Schwinger 
Equation (LSE) and discuss the obtained scattering am- 
plitudes and spectrum of the bound states. Finally, in 
Sec. HVI we generalize the discussion to the case of inter- 
acting many-systems including the cases of a ID degen- 
erate Fermi-g ID quasi-condensate and Tonks gas. 



II. MODEL 

In this section we introduce the our model system by 
specifying the Hamiltonian for a ID lattice gas coupled 
to an impurity, and we explain the key idea behind our 
setup. We will start with a discussion of spin-dependent 
collisions between the gas and the impurity, and then 
present the central idea of quantum interference as a way 
to switch atomic transport. 



A. Effective Spin-Dependent Hamiltonian 

We consider the dynamics of a spin- 1/2 atomic impu- 
rity Q coupled to a ID quantum gas of either bosonic or 
fcrmionic probe atoms A. The Hamiltonian for system is 
split into three parts as 



H = H A + H Q + H AQ . 



(3) 



Here Hq (Ha) describes the uncoupled dynamics of the 
impurity atom Q (the degenerate quantum gas of probe 
atoms A), while Haq accounts for the interaction be- 
tween the two atomic species, Q and A. 

A degenerate quantum gas of bosonic or fcrmionic 
atoms A trapped in the lowest band of a ID optical lat- 
tice is well described by a Hubbard model [6( 



H A = 



3 



(ij) 



U 
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where at (a,j) are the creation (annihilation) operators for 
an atom A on the site j, which obey standard commuta- 
tion (anticommutation) relations for the case of bosonic 
(fermionic) atoms A. Moreover, Eaj account for the 
shift of the bare energy of an atom localized on the site 
j in the presence of an external (e.g. magnetic) shallow 
trap, J is the tunneling matrix element for neighboring 
sites (ij) and U gives the collisional interaction, i.e. the 
onsite-shift for two atoms A localized within the same 
well (which would be zero for the case of fermions in 
the same internal state). Denoting the scattering-length 
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of the atoms A by a s , and their mass by m, we have 
U = Ai:h 2 a s J d 3 r\wj(r)\ 4 /m, where Wj(r) is the Wan- 
nier function for a particle localize on the site j. 

In the present setup we regard the impurity atom Q to 
be trapped within a tight one-dimensional lattice, as de- 
picted in Fig. ^b). Therefore, we may restrict ourselves 
to the lowest trap-state of the j = well for the internal 
states a =J.,T, respectively. The uncoupled dynamics of 
the impurity corresponds to spin-1/2 system, i.e. 



H Q = J2E Qta \Q a ){Q 



(5) 



where \Q a ) {Eq.v) denotes the state (energy) of the atom 
Q with spin a =|, |. 

Given the tight trapping of the impurity atom, the 
interaction of probe and impurity atom is restricted to 
the site of the impurity, an in general has the form of an 
effective spin-dependent collisional interaction 

Haq = W cSi1 \\)(Uala + W cSd \l}(l\alao, (6) 

where a\ (clq) is the creation (annihilation) operator 
for a probe atom on the site of the impurity, j = 0. 
Here, W cS ,a = 4irh 2 a a / ^ J d 3 r\w (r)\ 2 \^ Q ^(r)\ denotes 
the effective interaction for a probe atom A and the 
impurity atom Q in state a in terms of their effective 
scattering length a a and \i is the reduced mass for A 
and Q. The effective tunneling rate of a probe atom 
with energy E through the impurity is then given by 
J eff , cr = J /{E- WeB t(T ) +0(J 4 ) for the qubit in state a. 
An obvious way to provide for a spin-dependent single 
atom mirror is to have the effective interaction for one 
spin state as large as possible (| W e ff,T | J), thus block- 
ing the transport of the probe atoms through the impu- 
rity site, while for the other it is effectively not present, 
(|Wcff.j.| <§; J). This can be achieved, for example, by 
tuning the internal state dependent scattering length a CT 
or by engineering the spin-dependent trapping 6]. The 
quality of the qubit dependent switch then depends on 
the difference of the moduli of the effective interactions, 
IWeffjl — | Wcff.il- Thus the goal an efficient scheme is 
to make |W e ff -\\ — \ W e s i\ as large as possible and obtain 
|W eff , T |» J»|W efU |.' 



B. Controlling the transport by interference 

In this section we will show now that with the help 
of quantum interference we can engineer an effectively 
infinite (zero) atomic repulsion, W e gt,-f — > oo (W e ff,i — > 
0), for the qubit in state a —I (a =|). This is equivalent 
to tuning the Feshbach resonance to the point of infinite 
(zero) scattering length. 

The quantum interference mechanism required to engi- 
neer the described spin-dependence of W e fr, CT is obtained 
by exploiting the properties of either an optical or a mag- 
netic Feshbach resonance. In the case of an optical Fes- 
hbach resonance a Raman laser drives the transitions 




FIG. 3: The optical Feshbach setup couples the atomic state 
a \Q) (in a particular motional state quantized by the trap) 
to a molecular bound state \M) of the Born-Oppenheimer po- 
tential Vaq(t) with effective Rabi frequency il and detuning 
A. 



from the joint state of the two atoms on the impurity 
site, a\\Q a ) = \Q a ) £g> ajlvac), via an off-resonant excited 
molecular state back to a bound hetero-nuclear molec- 
ular state | Mo-) in the lowest electronic manifold (see 
Fig. |3J). The Raman processes is described by the ef- 
fective two-photon Rabi frequency fl a and detuning A CT 
for each spin-component a. For the case of a magnetic 
Feshbach resonance, the effective Hamiltonian has the 
same form, but with Q, a being the coupling strength be- 
tween the open and closed channels and A CT the detuning 
of the magnetic field. The Hamiltonian describing the in- 
teraction between the probe atoms and the impurity is 

m 

Haq = [E M ,<r\M a )(M a \ + fi CT (\M a )(Q a \a + h.c.) 

<7 

+W Q ,o|Qo)((9a|4«o + W A f,o|Mo)(M <7 |aJao] , (7) 

where the bare energy of the molecular bound state is 
Em, a = Ea,o + Eq,c + A CT . Here the first two terms 
describe the resonant coupling induced by the Fesh- 
bach mechanism, while the last two describe the off- 
resonant collisions between an atom A and an atom Q 
(a molecule M) in state a by means of their on-site 
shift Wq,o- (Wm,o) for the impurity site. The Hamil- 
tonian © conserves the spin-component of the impurity, 
S a = \Qa)(Q*\ + \M*)(M*\, i.e. [H,S a ] = 0. There- 
fore, we can consider the dynamics for the two spin com- 
ponents of Q separately, and in the following we will 
drop the spin index a and choose the reference energy 
as E A ,o = E Q , a = 0. 

For off-resonant laser driving (|A| ^> f2), the Fesh- 
bach resonance enhances the interaction between A and 
Q atoms, giving the familiar result W e g = Wq + tt 2 /A. 
However, for resonant driving (A = 0) the physical mech- 
anism changes, and the effective tunneling J e g of an atom 
A past the impurity (Fig. 21 I — > HI) is blocked by quan- 
tum interference. On the impurity site, laser driving 
mixes the states aJ|Q) and |M), forming two dressed 
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FIG. 4: A single atom passes the impurity (I — >III) via the two 
dressed states (II), |+) = oJ|Q) + \M) and |-) = aJ|Q) - \M) 
and quantum interference gives rise to an effective tunneling 
rate J e g. 



states with energies E± = Wq/2 ± (W|/4 + fi 2 ) 1 / 2 
(Fig. 01 II). The two resulting paths for a particle of en- 
ergy E destructively interfere so that for large Q, 3> J 
and Wq = 0, 



JcS 



J 



J 



E+n E-n 



o. 



This is analogous to the interference effect underl ying 
Electromagnetically Induced Transparency (EIT) |25| . 
and is equivalent to having an effective interaction 
W c ff — > oo. In addition, if we choose A = ti 1 /Wq, 
the paths interfere constructively, screening the back- 
ground interactions Wq to produce perfect transmission 
(W e ff — + 0). The insensitivity of the interference scheme 
to losses from the dressed states due to their large detun- 
ing has been argued in Ref. |18j . 



III. SINGLE PARTICLE SCATTERING FROM 
AN IMPURITY 

In this section we will analyze the scattering of a sin- 
gle probe atom A from an impurity atom Q. We will 
formulate the scattering problem, then solve the time- 
independent and time-dependent Schrodinger Equation, 
to finally obtain the dynamics of wave-packets in the lat- 
tice. 

We consider a probe atom A approaching the impu- 
rity from the left, as a plain Bloch-wave with quasi- 
momentum k. Hence the state of the system is given 
by 



(8) 



where \j) — \Q) ® aj|vac) is the joint state of the atoms 
A and Q, with A (Q) localized in the lowest vibrational 
state of the well j (the impurity well j = 0) and a is the 
lattice spacing. 



The free evolution of the system is given by the hop- 
ping of the atoms A between neighboring sites at the 
tunneling rate J, whereas the composite molecule M is 
detuned by A from the threshold for the joint state of 
A and Q. Thus, with Ek = — 2 J cos ka being the en- 
ergy of a Bloch-wave in the first Bloch-band with quasi- 
momentum fc, we have 

H ° = -J^(b' + l>01+h.c.) + A|M)(Af| 



7r / a 



7r / a 



dkE k \k)(k\ + A\M)(M\, 



(9) 



where \M) denotes the molecular bound state localized 
on the impurity site. From Eq. (O we obtain the propa- 
gation of the incoming plane wave \k) at group- velocity 
Vk — dEk/dk = 2Jasin£;a in the first Bloch band. 

Due to the strong confinement of the particles A and 
Q in the lattices, their interaction is restricted to the 
impurity site. There, their bare interaction induces an 
on-site-shift W for the joint atomic state of A and Q 
on the impurity (|0)). Moreover, the photo-association 
lasers effectively couple the latter state to the trapped 
molecular state (|M)) at Rabi- frequency O, yielding 



V = W\Q)(Q\ +fi(|M)(0| +h.c). 



A. Scattering solution 



(10) 



The scattering of a particle A with energy E = Ek in 
the first Bloch band by the impurity Q is described by 
a solution of the Lippmann-Schwinger Equation (LSE). 
The scattering wave function |</> + ) obeys 



- \k)+Go(E + tQ + )V\4> + }, 



(11) 



with incident plane wave | k) with quasimomentum k (0 < 
k < Tr/a), and Gq{z) = l/(z — Ho) the free propagator. 
Expanding the scattering wave function 



1/2 



5> 3 -|j)+/3|M) 



the amplitudes o>j and f3 satisfy 



, //••./ 



g 3 (E k )(Wa + np), 



(3 = g M (E k )ttao 
with atomic and molecular propagators 
9 3 {E) = (j\G o (E + i0+)\0) 



(12) 



(13a) 
(13b) 



a 

2^ 



+7r/a 



dk 



k | x j | 



Qm{E) 



la ^ - 

(M\G {E + i0+)\M) 



E — Ek + i0 + ivk/a' 



E - A + i0+ 
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Solving Eqs.lO we find 



a j 



Jk\xj\ 



P = 



iv k /a - W k 
—iQv k /a 



Q 2 + {E k - A + i0+) (W - iv k /a) 
with effective energy dependent interaction 

il 2 

W k = W- 



E k - A + i0+ : 



(15a) 
(15b) 

(16) 



where we read off the transmission and reflection ampli- 
tudes 



r k = 



1 



1 + iaW k /v k 
-1 

1 - iv k j aW k ' 



(17a) 
(17b) 



respectively. 

Note that the presence of the molecular state intro- 
duces an effective energy-dependent interaction W k . This 
can be interpreted in terms of an effective atomic scat- 
tering length with background scattering length propor- 
tional to W and a resonant term, corresponding to an 
optical Feshbach resonance at energy given by the detun- 
ing from the molecular state A, and width determined by 
the Rabi frequency f2. 

The scattering matrix 



S(E k ) 



r k t k 
t k r k 



(18) 



is unitary, as follows readily from the above expres- 



r k 



1. 



sions lfT7)). This implies T k + R k = \t k \ 
We can assign phase shifts t k ± r k = exp(i<S*) for 
the symmetric and antisymmetric states \k) ± | — fc), 
5~£ = — 2 arctan(aM / fe/wfc) and 5^ — 0, respectively, so 
that R k = sin 2 (£+/2). 



B. Discussion of the Scattering Amplitudes 

In the absence of molecular couplings (i! = 0) the on- 
site interaction W between the B atom and the impurity 
Q always gives rise to partial reflection and transmission, 
see Fig. |5ji, 



R, = 1 - T fc = 



W 



W 2 + 4 J 2 sin 2 (fca) 



< 1. 



(19) 



The significant new feature introduced by the optical 
Feshbach resonance is that we can achieve essentially 
complete blocking (R k = 1) and complete transmission 
(T k = 1). We obtain this in the limits fl^> J and A = 0, 
and 3> J and A = —£l 2 /W, respectively. Physically, 
the first case corresponds to tuning to the point of "infi- 
nite" scattering length, while the second case corresponds 
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FIG. 5: Reflection coefficient Rk, showing in (a) partial re- 
flection for bare background collisions W 7^ but Q = 0, i.e. 
qf = —1/2, —1, —2 (dotted, dashed, full line) and 7 = Er — 0, 
(b) complete reflection at fco = n/2a for 7/2J = 1/4,1,4 
(dotted, dashed, full line) and Er = qf = 0, cf. 7? ~ 1 for 
7/2J = 4 (full line), (c) perfect transmission at ko = n/2a 
for 7/2J = 1,4,16 (dotted, dashed, full line), qf — —2 and 
Er = 7/4, cf. R w for 7/2 J = 16 (full line), (d) the asym- 
metry of the Fano-profiles for qf = —1/4,-1,-4 (dotted, 
dashed, full line), 7/2 J = 1/2 and E R = 0. 



to tuning to the point of "zero" scattering length, respec- 
tively. 

In the general case the energy dependence of trans- 
mission and reflection has the form of a Fano-like profile 
(see Fig.[5j). In the region \E k \ <§C 2J we may neglect the 
dispersion effects, i.e. v k ~ 2Ja, and obtain Fano-line- 
shapes for transmission and reflection as 



T(e) 



R(e) = 



1 (e + q f y 



(e-l/qfY 



1 + 1 



(20a) 



(20b) 



where e = [E— Er) / '(7 /2) is the dimensionless energy in 
units of the resonance width 7, Er is the resonance en- 
ergy and qf is the Fano-q-parameter. These parameters 
of the Fano-profile are related to J, W, A, by 



Ajn 2 



7 = 
Er = 

If = 



W 2 
A - 

-W 
~2T' 



f 4J 2 ' 

wn 2 
w 2 + 4 j 2 : 



n 2 = J 7 (qj + 1) , (21a) 



A = E R -^, (21b) 



w 



-2Jq f . 



(21c) 



For W — the asymmetry parameter qj vanishes, the 
reflection profile is symmetric, and for £1 > resembles a 
Breit-Wigner-profile, see Fig[S^b). The maximum R k = 1 
is attained at e — —qf (E k = A, |A| < 2J) and has a 
width 7 = Q 2 /J. 
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For finite background collisions qf ^ (W 7^ 0) 
the transmission profile is asymmetric, and shows an 
additional minimum Rk — at e = —1/qf (Ek 
A - n 2 /W = A*, |A*| < 2J), see FigH^c). 

Near the edges of the Bloch band, k± — (tt =p 7r)/2a 
(£fe = ±2 J), transmission and reflection deviate from 
the Fano line-shape Eq. (|2(jp. There the group- velocity 
Vk — > and thus also the transmission Tk vanishes, unless 
the dressed resonance A* is tuned to respective edge of 
the Bloch band. The transmission coefficient are given 
by 



4 J 2 a 2 (k-k±) 



Ti 



k~k± 



[w - 



fl 2 /{AT2J)Y 
W 4 a 2 (k - k ± ) 2 
2714 



for A* ^ ±2J, 



for A* = ±2 J. 



The reflection coefficient Rk as a function of energy 
-2 J < E k < +2J is shown in Fig. El 

In the absence of molecular coupling, 7 = (O = 0), 
the reflection is unity at the band-edges, k — 0,n/a, and 
decreases within the Bloch-band due to the increase of 
the group- velocity (see FigEJa)). The profile is symmet- 
ric about the middle of the Bloch band, k = ir/2a, where 
it attains its minimum, 



Rk = 



1 



1 

-Qj 2 ~ 1 + {2J/W) 2 ' 



In the presence of molecular couplings, 7 (Q > 
0), and for En = qf = the reflection profile is still 
symmetric about k = ir/2a (see Fig^Jb)). However, now 
it approaches its maximum, Rk = 1, at k = it /2a, and 
has now two minima at k « 7r/4a and k « 37r/4a, given 
by 



i + (2j/ 7 ) 2 i + (V2j/n) 4 ' 



For 7, qj ^ we obtain an asymmetric Fano-profilc 
(see FigGJd)), which for \E R + q f j/2\ = |A| < 2J shows 
complete reflection, Rk = 1, at e = —qf, while for \Er — 
7/2g/| < 2J one has perfect transmission, Rk = 0, at 
e = +l/q f . 

The reflective and transmissive resonance are present 
regardless of the magnitude of qf, and their width is 
oc 7. However, for 7 < 8J|<7/ + l/<7/| they may both 
occur within the physical energy range of the Bloch- 
band (see FigOJd)), while for 7 > %J\qf + l/qf\ only 
one resonance appears (see Fig|SJb,c)). Thus in the 
limit 7 3> 8J\qf + l/<7/| we achieve complete block- 
ing , Rk = 1 for all k, by tuning En rj — (7/7/2 (see 
full- line in Fig|5fb)). Within the same limit we can 
also efficiently screen any background-interaction W and 
achieve complete transparency, X/. = 1 for all fc, by tun- 
ing £7r w l/2qf (see full-line in Fig.JSfc)). 



C. Interference mechanism 

Physically, the features of complete blocking (Rk = 0) 
and complete transmission (Tk — 0) are induced by an 
interference mechanism, as the probe atom may tunnel 
via two interfering paths of dressed atomic + molecular 
states, as depicted in Fig. 0] 

For simplicity we start by elucidating the underlying 
interference mechanism (present for f2 ^ 0) in the regime 
of strong coupling, O 2 3> (2 J) 2 + |WA|. In this regime 
we can consider the local dynamics within the individual 
sites, and treat the tunneling J by means of perturbation 
theory. 

To zeroth order in J the Hamiltonian H decouples the 
dynamics of the individual sites j as 

#(°) = W|0}<0| +ft|0)(M| + 

+ Q\M){0\ +A\M)(M\. (23) 

Outside the impurity (j 7^ 0) its eigenstates are the joint 
states of the atoms A and Q, \j), with energy Eq = 0, 
whereas on the impurity (j = 0) the strong coupling O 
between the atomic state |0) and the molecular state \M) 
induces the two states to split into two dressed state \E±) 
of atoms + molecules with energy E±, see Fig. QJ By 
diagonalizing the 2 x 2-matrix in Eq. (|23[) we obtain the 
amplitudes and energy of the dressed states as 



\E, 



E± = 



1/2 



W + A 



± 



|0)± 



n 2 



2 

( w - 



1/2 



\M), (24a) 



1/2 



W- A 



(W — A) + 4f2 2 



1/2 ; 



(24b) 



(24c) 



where £ characterizes the asymmetry of the amplitudes, 
i.e. for £ = (W — A) the dressing is completely 
symmetric while for |£| = I (f2 = 0) the atomic and 
molecular state decouple. From Eq. (|24b() we see that 
for Q, 3> |WA| the dressed states are far off-resonant 
from E = 0, and hence will be only virtually populated. 

The effects of the hopping of the atom A on the E sw 
modes \ j 7^ 0) can be accounted by means of an effective 
Hamiltonian i? c ff- Following Ref. pj we obtain the dy- 
namics as a perturbative series in the hopping amplitude 
J, H ca = + + + .... To first order in J 
one obtains 



H(1) = -j£|j>0'-i|- j£ljX-7 + 1 l + h - c - 



2 / dkE k \k a ){k a \ 



(25) 



a=L,R 



where |fci) (\kit)) are the Bloch- waves with quasi- 
momentum k on the left (right) side of the impurity site, 



L,R) 



/ a \ 1/2 
V27 



±j>o 



|J). 
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This the flat dispersion relation E^ = on the left and 
right side of the impurity is bent to e(k) — — 2J cos(ka), 
i.e. we recover the Bloch-band(s). 
To second order in J we obtain 



Jeff = "7T" 



-Jos 



J 2 A 



n 2 - wa' 



(26) 
(27) 



We see that tuning on resonance A = the two contri- 
butions in Eq. (|27|l cancel each other as 



Jeff = 



J 



P 



Vw 2 + n 2 vw 2 



n 2 



= o, 



(28) 

which gives perfect blocking by the impurity Further- 
more, from Eq. I|24bjl we obtain that for A* = A — 
il 2 /W = one of the dressed states \Eb) becomes a 
resonance for an incoming particle {Eb = Eq = 0) and 
for £1^0 provides for complete transmission by means 
of photo-assisted tunneling. The described interference 
mechanism induced by the optical Feshbach resonance 
is in marked contrast to the situation where one has 
background collisions. There the particle A can tunnel 
only via one path through the impurity (|£| = 1), and 
therefore the effective hopping rate is always finite, i.e. 

Jeff = ~J 2 /W ? 0. 



D. Discussion of Bound-states 



T(z) 




® 


pole 




cut P° le 


E_ 


-2J 


+2 J E + 




W/2J 




(d)+5 



E/2J 




FIG. 6: (a) The analytic properties of the T-matrix in the 
complex plane. The cut on the real axis is due to the prop- 
agation of (dressed) Bloch-waves. The two poles on the real 
axis correspond to two bound states and give rise to interfer- 
ence, (b-d) The spectrum E n of the Hamiltonian H = Ho + V 
as a function of (a) the on-site-shift W on the impurity for 
fi = 0, (b) the detuning A for ft/ J = 2 and W = (solid 
lines), for fi/J = 4 and W = (dashed lines), for Q/J = 2 
and W/J — 4 (dash-dotted lines), (c) the Rabi-frequency Q 
for A = W = (solid lines), for A/ J = 4 and W = (dashed 
lines), for A/ J = W/J = 4 (dash-dotted lines). The con- 
tinuous (scattering) spectrum, i.e. the Bloch-band Ek, is in- 
dicated as a shaded region — 2J < Ek < 2J, whereas the 
bound-state(s) Eb are shown as (solid, dashed, dash-dotted) 
lines. 



For completeness we here derive the exact bound-state 
spectrum of H. For the exact scattering solution, de- 
tailed in Sec. IIIII the bound-states take the role of dressed 
states \E±), which are responsible for the interference 
mechanism. We will show that for arbitrary £1 there are 
always two bound-states, provided |A — 2 /IT| > 2J. 
For | A — ri 2 /TT| < 2 J one of the bound-states turns into 
a resonance, which makes the impurity completely trans- 
parent for the atom A, Tfc = 1. Furthermore, we will 
show that the bound-states for extends over several lat- 
tice sites for Sl,W,A ~ J. This is in marked contrast 
to the perturbative result, where the dressed states were 
localized on the impurity site, cf. Eq. I)24a[l . 

We obtain the bound states wavefunctions \4>b) from 
the homogeneous Lippmann-Schwinger equation 



,) = G {E)V\<j> B ), 



(29) 



where <pB denotes the bound-state with energy E = Eb 
(\Eb\ > 2J). Using the ansatz 



\4>b) = 5>ili>+/W 



(30) 



we find that the atomic and molecular amplitudes, ctj 
and f3, satisfy 

<h = Gj{E B )[Wa +£ll3], (31a) 
P = g M (E B )^ao. (31b) 



The atomic and molecular propagators are given by 

e-M/rs [ S ign(-£; s )] j 



Qi{E B ) = (j\G Q (E B )\0) = 



-2 Jsinh (kbo) 



Qm{E. 



B 



(M\G o (E B )\0) 



1 



Er 



and r_B = a/acosh(|i?B|/2J) > denotes the size of the 
bound-state. 

For convenience we first consider the case O = 0. There 
the molecular state decouples from the atomic ones, and 
we have one bound-state \<j)\) = \M) with energy Ei = A. 
Moreover, for W ^ we have another bound-state \4>2) 
with energy E 2 = sign(W)^W 2 + (2 J) 2 . Its amplitudes 
are given by (3 = and 



tanh ( — 

tb 



"I 1/2 



o-\xj\/r B 



sign(-VF)] J 



with the size rs = a/arccosA/l + (IT/2J) 2 . The spec- 
trum of the system is plotted in Fig. EOT) as a function 
W . For attractive (repulsive) interaction W < (W > 0) 
the energy of the bound-state </>2, lies below (above) the 
Bloch-band, i.e. E 2 < — 2 J (E 2 > 2 J), respectively. For 
\W\ < 2 J bound state 4> 2 extends over several lattice 
sites, while for \W\ ^> 2 J it is localized on the impurity. 
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In the case f2 7^ a nontrivial solution of Eq. 1311) 
requires 



23 



1/2 



W 



n 2 



Eg -A' 



(33) 



which determines the bound-state spectrum 
Eb(W,A,£1). From Eq. IpTT)) we obtain the atomic 
and molecular amplitudes as 

^ ^B-A e _ |g . |/rB [gign( _ jBs)] , ^ (34a) 



SI 



/3 = 



-1/2 



(34b) 



We solve Eq. Ij33(l by expressing one of the parameters, 
either W, SI or A, in terms of the bound-state energy 
E B , 



A(E B ) = E L 



n 2 



E Bx /l-4:J 2 /E 2 B + W' 



&{E B ) 



a 2 



/ 4 7 2 

W ^ =E °f + E B - ' Eb-A- 



(35a) 

1/2 
(35b) 

(35c) 



Inverting the functions Eq. I|35|l for Eb yields E B as a 
function of A, fl and W, respectively. For fixed A, Q, 
W we carry out the inversion by plotting in Fig. the 
r.h.s. ofEq. I|35|) ( a.b.cl as a function of E B ■ In particular 
in Fig. E{b) we plot detuning A as a function of Eb for 
constant f2 and W, in Fig. EJc) we plot Rabi-frequency 
Omega as a function of Eb for constant A and W, and 
in Fig. HJd) we plot on-site shift W as a function of Eb 
for constant A and SI. In the following we will give a 
detailed discussion of Fig. 

For no background collisions, W — 0, and arbitrary de- 
tuning A, one always has two bound-states 01,2 with en- 
ergy Ei < — 2J and E2 > 2J, respectively, see solid line 
in Fig. EJc) corresponding to SI = 4 J. For |A| >• fi the 
energy one bound-state approaches A and is wavefunc- 
tion becomes localized on the impurity, while the energy 
the other approaches the Bloch-band and its wavefunc- 
tion extends over several lattice-sites. For A = the two 
bound-states are split symmetrically, and their energies 
are given by (see solid line in Fig.HJd)) 



E h2 = ±\/2J 2 + V4J2 + fi2. 



(36) 



The symmetric splitting allows for complete reflection at 
Ek = 0. In the limit SI ^> 2 J we recover the perturbative 
result Eq. I|24b|) . as the energies of the bound states ap- 
proach E± = ±S1, with their wavefunctions given by the 
dressed states \E±), cf. Eq. I|24a|l . 



For finite onsite shift, W 7^ 0, we also have two bound- 
states provided |A*| = |A— S1 2 /VF| > 2 J, see dashed lines 
in Fig. EIc,d). With increasing detuning A the energy 
of the bound-state tfii approaches the Bloch-band from 
below until crossing it for -2 J + ft 2 /W < A < +2J + 
tt 2 /W. In this parameter regime there is merely one 
bound state, while the other develops a resonance. This 
allows for perfect transparency at Ek = A — SI 2 /W . 

Finally we remark that the size rs of the bound-states 
is inversely proportional to the separation of their en- 
ergy Eb from their Bloch-band. Thus for \Eb\ -C 2 J 
the wave-function of the bound-states extends over sev- 
eral lattice sites. For \Eb\ 2 J the bound-states are 
localized on the impurity and we recover the results of 
the previous section. 



E. Wave-packet dynamics 

As an illustration of the time-dependence of the in- 
terference mechanism we simulate the evolution of a 
gaussian wave-packet ip(t) with mean quasi- momentum 
k = 7r /2a incident from the left of the impurity. 

These wave-packets are obtained as superposition of 
the scattering solutions </>+ of Sec. IIIII i.e. their atomic 
and molecular amplitudes, Oj(t) = (j\ip(t)} and (3{t) = 
(M\ip(t)), are obtained as 

= (37a) 
3' 

f3(t) = J2u Mtj > (*W(0), (37b) 
3' 

with the full propagator for the system given by 

e -(\Xj>\ + \xA)/ r B-iEBt s ig n (_E B y'+j 

B 



U 



3', 3 



(*) = £■ 



2tt 



dke~ iE 



\ e +ik\x } t-x } \ _|_ rke + ik (\ x j' \+\ x i\) 



e -M/^- iE Bt s i g n{-E B y , 

UM,j(t) = = h 



1 , \E B \(E B -A)2 



/ dke~ iEkt 

2tt 



nt k 



+ik(\ Xj \) 



Ei. 



i0^ 



(38) 



On the left side of Fig. [7| we plot the atomic popu- 
lations of the individual sites, n(xj,t) — \aj(t)\ 2 . The 
right side shows the corresponding atomic populations of 
the atom A on the left, n L = Ylj<o n(xj,t) (dashed line), 
and on the right side of the impurity, tir = J2j>o n(xj,t) 
(dashed-dotted line). We also plot the population on the 
impurity, i.e. the atomic population, hq = n(xo,t) (solid 
line), and the molecular population, njvf(i) = K-^IV'W)! 2 
(dotted line) . The three different sets in Fig. [7| corre- 
spond to different coupling strengths, SI, and detunings, 
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FIG. 7: Evolution of a Gaussian wavepacket with mean quasi- 
momcntum k = n/2a for W — 2J and (a) £1 = 0, (b) = 2J, 
A = 0, and (c) Q — A — 2 J. On the left side the atomic den- 
sity, n(xj, i), of atoms A is plotted (darker region correspond 
to higher density). The right side shows the corresponding 
population of atoms A of the sites to the left, nz,(i) (dashed- 
line), to the right, nn(t) (dashed-dotted line), and on the 
impurity site no{t) (solid line), as well as the population of 
the molecular state, nu{t) (dotted line). 



A. For all cases we choose W — 2 J. In Fig. Ufa) we have 
£1 = 0: the atom is partially reflected from the impurity 
with R k = 1. In Fig. □» we set = 2J and A = 0, 
which gives rise to complete reflection of the wavepacket, 
Rk = 1. In Fig. EKc) we have fl = 2 J, but now A = 2 J. 
We have complete transmission of the atom through the 
impurity, = 1. All this is consistent with the results 

of sec, inn 



IV. MANY BODY SCATTERING FROM AN 
IMPURITY 

In this section we will analyze the evolution of a ID 
lattice gas of many atoms A interacting with an impu- 
rity atom Q. Since the statistics of the atoms A plays 
a dominant role, we will consider the cases of fermionic 
and bosonic atoms, separately. In this context we will 
study analytically the limiting cases of an ideal Fermi- 
gas, an ideal Bose-gas and a Tonks-gas. An exact numer- 
ical treatment of the dynamics for the lattice-gas having 
arbitrary interaction U is given in Ref. |24j . 



A. Ideal Fermi-gas 

We first consider the case, where the probe atoms A are 
spin-polarized fermions. The Hamiltonian for the system 
is given by 



% ]+l a 3 



A\M){M\ 



+Q[\M)(Q\a + al\Q)(M\^ 
+W Q \Q){Q\ala + W M \M)(M\ala a , 



(39) 



where the operators Oj (ctj) create (annihilate) an atom 
A on site j, and obey the canonical anti-commutation re- 
lations {a;, ctj} = Sij and {a*, a 3 - } = {ctj, at} = 0. More- 
over, \Q) (\M)) denote the states with an atom Q (a 
molecule in state M) on the impurity, and Wq (Wm) is 
the onsite-shift for an atom A and an atom Q (a molecule 
M) on the impurity. 

For simplicity henceforth we will restrict ourselves to 
the case of equal on-site shifts Wm = Wq = W. In this 
case we may rewrite the Hamiltonian as 

H = — faja, + i + aj +1 ajj + Wala + 
j 

A/t/ + n(/t ao + 4/), (40) 

where the ladder operators f> = \M)(Q\ and / = \Q)(M\ 
obey standard fermionic anti-commutation relations and 
anti-commute with aj and aj. The corresponding equa- 
tions of motions for a,j and / are linear, provided Wm = 
Wq. Thus for a Fermi-gas of N atoms A the scatter- 
ing off the impurity atom Q will occur independently for 
each fermion afe with scattering amplitudes t& and r^, 
according to their quasi-momentum k, cf. Eq. (fTTjl . The 
details of this calculation will be given below. 

We will here detail the time-dependent scattering for 
a Fermi-gas of N atoms A. For concreteness we assume 
the fermions to be initially trapped in a box of M sites to 
the left the impurity atom Q. The corresponding wave- 
function of the system is given by 



N 



i*(*=o)>=n 



\Q), (41) 



Otjikn) 



FY r sm(k n Xj) for - M<j < 1, 







else, 



where the quasi-momenta k n = nir/(M + l)a. This cor- 
responds to a Fermi-sea filled up to Ep = 4J sin 2 (k Fa), 
where kp = vn /a is the Fermi-momentum and the initial 
filling-factor v = N/(M + 1). 

At time t = we open the impurity (cf. Fig. ^b)), and 
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FIG. 8: Evolution of a Fermi-gas with filling factor v = 1/2 
(N = 25 particles on M = 50 sites) after opening the switch 
at t = 0. The effective Rabi-frequency is (a) Q = 0, (b) S7 = J 
(c) f2 = 2J, and in all three cases we have W — A = 0. The 
left side shows the atomic density n(x, t) (darker regions cor- 
respond to higher density), and the right side the momentum 
distribution n(k,t) at time t = 0,M/2J, M/J (from bottom 
to top). The corresponding density distribution at this three 
times are indicated by arrow. 



from Eq. H40|) we obtain the evolution of the system as 



N 



!*(*)) = n 



\Q), (42a) 
a ] (k n ,t)=Y,U^(t)a jl (k n ), (42b) 
I3(k n ,t) = U M j> {t)a r (fc„), (42c) 



where U a .j(t) are the single-particle propagators, cf. 
Eq. . According to Eq. (|42a|l the scattering from the 
impurity occurs independently for each particle in the ini- 
tial Fermi sea, with scattering amplitudes tk and given 
in Eq. (|17fl for < k < kp. The atomic and molecular 
densities are thus given by the sum of the probabilities 
for the single fermions in the Fermi-gas, 

N 

n(xj,t) = {a) j a j )t = S ^\oi j {k n ,t)\ 2 , (43a) 

n=l 

N 

n M (t) = l-n Q (t) = Y,\(3(k n ,t)\ 2 . (43b) 



Moreover, after opening the switch the atomic quasi- 
momentum distribution in the Fermi-gas for the semi- 
infinite system is given by 



n(k, t) 



2tt ^ 



(at, a,) 



JV 

a x -> 
2^^ 

n=l 



j 



'Ctj(k n ,t) 



(44) 



In Fig. IHa,b,c) we show the evolution for a Fermi-sea 
with v fv 3/4, i.e. N — 38 particles initially on M = 50 
sites. For each simulation we have W — A = 0, but 
the driving varies as fl/J — 0,1,2 in Fig. |S[a,b,c) re- 
spectively. On the left side we plot the atomic density 
n(xj,t) (darker regions correspond to higher density) . To 
the right we plot the respective momentum profiles n(fc, t) 
for the Fermi-gas at times t = 0,M/2J, M/J (from bot- 
tom to top). This times are indicated by arrows each 
figure. 

In Fig. IHIa) we see the evolution of the noninter act- 
ing system, O = 0. The atomic cloud expands freely to 
the right after opening the switch at t = 0. The cor- 
responding momentum distribution is initially given by 
n(k,t = 0) w 9{k F - \k\)Ma/2ir (see profile at the bot- 
tom). With progressing time < t < M/2J the gas 
develops a forward peak at k = it /2a (see profile in the 
middle) until becoming asymmetric as n(k 7 t = M/J) m 
(M/2n)8(k F - k)29{k) for t < M/J (see profile at the 
top). 

In Fig. Efb) we show the behavior for for weak laser 
driving, f2 = J. We notice that there is already substan- 
tial blocking by the impurity. The corresponding mo- 
mentum profiles show that the blocking is mainly due to 
the complete reflection of fermions with quasi- momentum 
near k = ir/2a to k' = —it /2a. 

In Fig. |5f c) we plot the densities for resonant driving 
with £7 = 2J. The transport through the impurity is 
efficiently blocked by the impurity atom, as the initial 
densities n(xj,t = 0) and n(k,t = 0) are almost com- 
pletely preserved. 

In the following we will consider the number of parti- 
cles on the right side of the impurity, 



N R {t) 



j>0 



n(xj, t), 



(45) 
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FIG. 9: The number of atoms A transmitted through the im- 
purity as a function of time t, Nn(t), for a Fermi-gas with 
v = 1. (a) We gave A = W — and the coupling is Q — 
(solid line), fl = J (dashed-dotted line), f2 = 2 J (dashed 
line). In (b) there is no coupling, fl — 0, but the background 
interaction are W = 0, J, 2 J, 4J (solid, dash-dotted, dashed, 
dotted line). The system establishes a constant particle cur- 
rent up to t ~ 50/ J A = W = and the coupling is = 
(solid line), Q — M/J = 50/J, and saturates until all N = 50 
atoms A passed the impurity. 



and the corresponding particle current through the im- 
purity, Ir. 



(46) 



They characterize the behavior of the switch. 

In Fig. |5Ja) we show the number of particles Nji(t), 
for the same parameters as in Fig. |H1 i.e. for each the 
initial filling factor v w 3/4 and W = A = 0. The solid 
line shows Nr for the no coupling to the impurity, f2 = 

0, and corresponds to the densities shown in Fig. Ufa). 
The dashed-dotted line shows the behavior for Q = J, 
cf. Fig|Sfb), and the dashed line corresponds to fl — 2 J, 

1. e. FigOJc). Moreover, in Fig.[5Jb) we show the number 
of particles N_n(t), for initial filling factor v « 3/4 and 
il = A = 0, but for an onsite shift W = 0, J, 2 J, 4 J (see 
solid, dash-dotted, dashed, dotted line), respectively. In 
general, after a short transient period, of the order of the 
inverse tunneling rate 1/J, the number of particles on 
the right side of the impurity increases linearly with t. 
Thereby the system establishes a roughly constant flux 
of particles Ir through the impurity. The flux persists 
up to t w M/J, which is indicated by a vertical dotted 
line in Fig.|5fa,b). Then the population on the left side 
of the impurity is substantially depleted and therefore 
Nji(t) saturates until all particles tunneled through the 
impurity, yielding NR(t) = N and Ifl(i) = for t — > oo. 

We are interested in the linear regime. From Eq. I|43a|) 
we obtain the constant average current as (cf. the 
Landauer- formula) 



dN h 
dt 



1 

27 



dkv k T k , 



(47) 



where Vk — — 2Jasin(fca) is the group-velocity of the 
quasi-particles with quasimomentum k, kp — vn/a — 
Nit/{M + \)a the Fermi-quasimomentum and Tk the cor- 
responding transmission coefficients, cf. Eq. i|17fl . Thus 




FIG. 10: The quasi-steady state atomic current In through 
the impurity. Its dependence on the filling factor v is shown in 
(a) for Q = but W = 0, J, 2J, 4J, and in (b) for W = A = 
but O = 0, J, 2 J. The dependence of Ir on the detuning A is 
shown in (c) for Q = 2J, W = and in (d) for fi = W = 23. 



the average current is obtained by integrating the Fano- 
profiles Tfc = 1 — Rf., see e.g. Fig up to the Fermi- 
momentum. 

For an uncoupled impurity, W = il = 0, we have Tk = 
1. Thus the current is given up to a constant by the 
Fermi-energy, 



r(o) 



27 



2J 



U , n / VK \ 

— sin — 

7T V 2 J 



(48) 



In Fig. IHtb) we plot the dependence of iS' as a function 
of the filling factor v as a solid line. 

For a finite on-site shift, > 0, but no laser-driving, 
= 0, we have T k < 1, see Eq. JT^J and Fig. EI a). Thus 
the current through the switch decreases as 



AW) 



El 

2tt 



w 2 , fw 2 

arccoth 

2ttE V E E F 



2JEF ) (49) 



where E Q — y/W 2 + 4 J 2 is the modulus of the bound- 
state energy. The exponential decay of Ir^ with in- 
creasing coupling strength W, i.e. the arccoth-term in 
Eq. 1(35} . is characteristic for a system with one bound- 
state. The dependence of the current I R W ^ on the filling- 
factor v is shown in Fig. HOf a). The solid line shows 
the non-interacting value, W = 0, while the dash-dotted, 
dashed, dotted line correspond to W = J, 2 J, 4 J, respec- 
tively. 

However, for resonant driving at > and couplings 
W = A = we obtain a symmetric Fano-profile for Rk 
with respect to k = n/2a, see Fig. HJb). Therefore, by 
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integrating the latter profiles we obtain the current as 



for the system is given by 



r (n) 



Ep E+E± 
2^ ~ 2tt (El - E 2 _ 

ElE- 



-arccoth 



E 2 _ + 2JE F \ 



Ep 



2tt Ei 



-arccot 



E 2 +2JE F \ 



E— Ep 



) 



(50) 



with E± = VV^ 4 + 4J4±2J 2 . The arccoth-term in 
Eq. I|5U|) gives the mean effect of the reflection as that 
typical of a system with one bound-state, while the os- 
cillating arccot-term is induced by the presence of two 
interfering poles in the scattering matrix. The current 
Eq. (|5U|I is plotted in Fig. llUf b') as a function of the initial 
filling factor v for the uncoupled impurity, O = (solid 
line), for fl = J (dashed-dotted line), and for = 2 J 

(dashed line) . For finite driving the current 1^ shows a 
plateau at v — 1/2, as all the Bloch-waves near k = it /2a 
are completely reflected from the impurity (see Fig.[5jb)). 
From Eq. I|50|l we obtain that already for f2 ~ 4J the 
current of particles through the impurity is completely 
suppressed for arbitrary filling v, i.e. up to Fermi-energy 
Ep = 4 J. 

In the following we discuss the dependence of the cur- 
rent for Vi > on the detuning A. In Fig. HOf c) we 
show the current Ir for £1 = 2 J but still W — as a 
function of the detuning A for several initial densities v. 
The solid line corresponds to commensurate initial filling, 
v = 1, the dash-dotted line to half-filling, v — 1/2, and 
the dotted line to a dilute Fermi-gas with v = 1/4. The 
current shows a symmetric profile with a minimum at 
A sa — 2Jcos 3 (^7r/2) and approaches its threshold value 

IjP for | A | 3> f2. Notice that the resonance for the many- 
body Fermi-gas with increasing density from the bottom 
of the Bloch-band, A w —2 J, toward the middle of the 
band, A = 0. 

For finite W = 2 J (see Fig. HOf d)) the dependence 
shows an asymmetric profile and reaches its threshold 
value 1^ (cf. Eq. gSJ) for large detuning, |A| > 
2 /| W|. We notice that although the single fcrmions 
in the Fermi-sea scatter independently, we obtain a finite 
current for Q < J, even on resonance. This is caused by 
the fact that the various fermionic modes see the res- 
onance a different (energy-dependent) detuning A — Ek, 
which leads to a shift of the minimum (and maximum) 
of the transmitted current proportional to the density of 
the Fermi-gas, see Fig. HUf c.d'). However, in the limit 
of strong driving f2 3> J we recover the features of per- 
fect blocking for A ss and of perfect transmission for 

a w n 2 /w. 



B. Ideal Bose-gas 

We now consider the case, where the probe atoms A 
are spin-less non-interacting bosons. The Hamiltonian 



H = + h - c -) + A|M)(M| - 

j 

+W Q \Q)(Q\ala + W M \M)(M\ala 
+n(\M)(Q\a + al\Q)(M\) , 



(51) 



where the operators (a,j ) create (annihilate) an atom A 
on site j, and obey the canonical commutation relations 
[aj,at] = Sij and [a^et.,] = [aj,aj] = 0. Moreover, \Q) 
(|M)) denote the states with an atom Q (a molecule in 
state M) on the impurity, and Wq (Wm) is the onsite- 
shift for an atom A and an atom Q (a molecule M) on the 
impurity. As in the previous section we will henceforth 
restrict ourselves to the case of equal on-site shifts Wm — 
Wq = W . In this case we may rewrite the Hamiltonian 
as 



H 



—J f a;taj-+i + a] +1 aj ) + Wa^ao + 

3 

+Aa + <J~ + I er + ao + a a~ 



(52) 



where the pauli operators er + = \M)(Q\ and a = 
\Q)(M\ obey canonical anti-commutation relations and 
commute with a,j and at. The Hamiltonian l|52l) corre- 
sponds a multi-mode Jaynes-Cummings model. 

In the following we will consider the scattering of a 
gaussian wavepacket of iV bosons A, all initially occupy- 
ing the same single particle state, a(xj, t — 0), approach- 
ing the impurity atom Q with mean quasi-momentum 
it ja > fco > and width <5fco <C it /a. The corresponding 
wavefunction for the system is given by 



|*(t = 0)) = |Q>- 



1 



a(xj, t = 0)aj 



iV 



|vac)[53) 
(54) 



where for 5ko -C it /a the normalization is given by Af 2 = 
{25k^a 2 /it) 1 / 2 = uq/N in terms of the peak density of the 
gaussian wavepacket, n = n(x(0),t = 0), and x(0) <C 
denotes the mean position of the particles A at t = 0. 

For Q — the equations of motion for aj decouple from 
ct~ . Therefore, we obtain the scattering of the bosons A 
by the impurity, as 



!*(*)) = 10)- 



N 



(55) 



where the single-particle wavefunction for finite W was 
already obtained in Sec. IIIII For this case all the re- 
sults obtained in Sec. IIIII hold, and we obtain e.g. the 
density as N times the single particle result, n(xj,t) = 
N\a( Xj ,m 
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FIG. 11: Scattering of a Gaussian wavepacket of bosonic 
atoms B with mean quasi-momentum ko and width Sko as 
obtained by linearizing the impurity (see text), (a) Molecu- 
lar population as a function of time for a narrow wavepacket 
with k = tt I a and Sk = 7r/100a for ft = J, J/2, J/4 
(solid, dash-dotted, dashed line), (a) The maximal attained 
molecular population , n* m , during the interaction process 
with a weakly coupled impurity (fi = J, A = W = 0), 
as a function of the quasi-momentum ko of the incoming 
wavepacket. The three lines correspond a wavepacket with 
width Sko = n /200a, n /50a, iv /20a (solid, dash-dotted and 
dashed line). The dependence of n* m on the parameters Q 
and A for k = n/2a and Sko = vr/50a in (c) without the 
presence of background collisions, W = 0, and in (d) for an 
onsite-shift W = 2J. 



1. Linearization of the impurity 

For fl > J > |A|, |W|, we obtained in Sec. Ill 
that the population of the molecular state was strongly 
suppressed, i.e. as (J/il) 2 , and thus we approximate 
a z — > —1. Thus we linearize the spin, i.e. set a + — > 6^ 
and <7 - — > b, where b and &t obey canonical commutation 
relations. The scattering of the bosons A by the impurity 
is given by 



!*(*)} 



i 



N 



|vec)a,(56) 



where the single-particle wavefunction ctj (t) and the am- 
plitude (3{t) of the molecular state 6^|vac) = <j + \Q) — 
\M), were obtained in Sec. HTT1 (see Eq. 

Self-consistency of the replacement a z — > — 1 requires 
that the obtained molecular population n rn (t) <C 1. From 
the linearization we obtain the molecular population as 
(see App. EJ 



M = 



dk — e 

E k - A 



fc 2 /45fc^-ifcs(0) 



(57) 



analytically e.g. by using a saddle-point method, see 
App. [X] We find that the maximal attained molecular 
population, n* n , is proportional to the initial density no 
of the gas. In the case of a broad resonance, f2 > Sko, we 
find 



n Q v 2 /a 



(E - A) 2 v 2 /a 2 + n 2 + W(E 



Ay 



r(58) 



with Eq = — 2Jcos(fco&) and f o = 2Jsin(/coa)- Moreover, 
for a extremely narrow resonance, £1 < 5ko and |A| < J, 
one obtains 



n 



n 2 



25k 2 a W 2 + v 2 /a 2 



-(ka-k„) 2 /28k„ 



(59) 



where k = k — ko, and the Fourier-integral is obtained 



where v+ — 2Jsin(/c*a) and k+ w arccos(— A/2 J) denotes 
the position of the maximum of Tk/(Ek — A) 2 . 

In Fig Jllf a) we show the molecular population, n m , as 
obtained by the replacement cr z — > — 1. In Fig. Illf a) we 
plot the n m (t)/no for an incoming gaussian wavepacket 
with ko — it /2a and 6k Q = 0.0l7r/a for driving O = J/4 
(dashed line), Q = J/2 (dashed-dotted line) and CI = J 
(solid line). In all three cases we have A = W = 0. 
We see that with increasing Rabi-frequency fl the at- 
tained molecular-population quickly drop as ~ J 2 /il 2 , 
and that the molecular-population closely resembles the 
density distribution n(xj,t = 0) of the atomic-cloud A. 
In Fig. Illf bl we plot the maximal attained population, 
, as a function of the incoming momentum of the gas, 
ko, for fl = J and W = A = 0. The four lines cor- 
respond to different width Sko of the wavepacket, i.e. 
Sko = 0.00l7r/a (solid line), Sko = 0.0l7r/a (dashed 
line), Sko — OAn/a (dash-dotted line) and Sko — 0.2ir/a 
(dashed line). For a given width Sko the molecular 
population attains its maximum for fco = n/2a, i.e. 
where we have complete reflection of the wave-packet (see 
Fig-EIb)). At the point of complete-reflection, fco = 7r/2a 
the population attains its overall maximum for a nar- 
row momentum-distribution, i.e. for Sko ~ * we have 
n m « 4n for O = J and A = W = 0. The depen- 
dence of n* n on the detuning A and the Rabi-frequency 
fl is shown in Fig. UTT cd) for W = and W = —2 J, re- 
spectively. In both figures the gaussian wavepacket has 
fco = n/2a and Sko = 0.027r/a, i.e. initially extends about 
~ ir/aSko = 50 lattice sites. For W — we have complete 
reflection of the wavepacket for A = 0, and the attained 
molecular population is maximal, n m w 3noi for Q w 2 J 
and A«0 (see Fig. lTTT c)). However, for finite W = —2J 
we have also complete transmission of the wavepacket, 
i.e. for A = Q./W. From Fig. Illf d) we notice that for 
a given Omega the maximal population is shifted from 
A«0 towards the point, where one has complete trans- 
mission of the wavepacket, i.e. A rj Q/W. However, 
the overall maximum of n* n in both cases, W = and 
W = —2 J, is attained for f2 « 2 J and for stronger driving 
quickly drops as n* m » no(2J/fl) 2 . As the replacement 
<7~ — > & is self-consistent for a dilute gas with densities 
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no <C (nj^/no) -1 , we see from ^2 that the approxima- 
tion holds, even on resonance A = 0, for strong driving 
57 > 4J up to densities as high as no ~ 5 and for small 
densities no < 1 only fails for A ss and f2 » 2 J. 



2. Time- dependent Variational Ansatz 

In the following we use a time-dependent variational 
Ansatz to describe the behavior of the many-body wave- 
function in near resonance A sa for ~ J, i.e. in the 
regime where the approximation a z — > — 1 fails already 
for small densities, no > 0.05. As a generalization of 
Eq. (|30(l for N > 1 bosonic atoms A we choose as an 
number-conserving Ansatz for the state of the system 



|*(i)> = CQ(t)\Q) 



M*) 1 



i N 



c M {t)\M) 



vw. 

[a M (t) r 



vac) 



N-l 



where a a (t) — J2j a j,<r(t) a j represent two 
orthogonal time-dependent modes for the field of the 
bosonic atoms A given that the impurity is in state 
(7 = Q,M, and the amplitudes for the impurity, c CT (t), 
and for the bosonic wavepackets, otj,<r(t), are normalized 
as 12a | Co- (t) | 2 = J2j \ a j,<?(t)\ 2 - The equation of motion 
for variational parameters, c a {t) and ay l<7 (t), are obtained 
by minimizing the corresponding action (cf. App. [Bj 



(60) 




FIG. 12: Evolution of a Gaussian wavepacket with mean 
quasi-momentum k — n/2a and width Sko — 0.02-7r/a ob- 
tained from the variational Ansatz Eq. I6UI . In (a) the re- 
flection coefficient R is plotted as a function of the Rabi- 
frequency fl for A = W = 0, and (b) shows the corresponding 
attained maximal molecular population n* m during the inter- 
action with the impurity atom. The reflection coefficient as a 
function of the detuning A for Q. — J and (c) W = and (d) 
W = 2 J. The various solid lines correspond to initial densities 
n = 0.1, 0.15, 0.20, 0.25. The dotted-line is for n = 0.05 and 
corresponds to a linearization of the impurity, see Eq. 1561 . 



S(t) = 



(tt t |W t )-QP t |W t 



(*t\H\V t 



(61) 



with respect to c„(t)* and as a set of coupled 

non-linear differential equations, which we integrate nu- 
merically. Thus we obtain the dynamics of the system. 
In Fig. E| we show the obtained reflection-coefficient 



R = lim — n(xj , t) 
t-koo N v ■" ' 

j<0 



(62) 



and the attained peak molecular population for 
a Gaussian-wavepacket with narrow momentum fco — 
tv /2a, i.e. Sko — 0.02n/a <C 7r/a for initial density 
n = 0.05,0.1,0.15,0.20,0.25 (values indicated in plots). 
In Fig. H2» (Fig. 0b)) we plot R ( function 
of the Rabi frequency for A = W = 0, i.e. when 
complete reflection of the wavepacket was predicted by 
the bosononic approximation of the spin. The reflection 
shows a non-linear behavior in the density for > .1, i.e. 
decreases as ~ 1 /no with increasing density no. While for 
no = 0.05 (see dotted line) we have complete reflection of 
the wavepacket for fl > 2 J, for higher densities we still 
have a finite transmission at £1 = 2 J. From Fig. I12f a1 
we see that with increasing density the transmission co- 
efficient rapidly deviates from the low(zero)-density re- 
sult and approaches a linear behavior in Q already for 



no = 1/4. In Fig. ll2f a) we plot the dependence of n* m on 
O, which shows that the maximal population is attained 
for w J/2 and decreases as ~ (J/O) 2 for fl ^> Jno- 
Moreover, we see that with increasing density, no > 0.10, 
the peak in the molecular population, n^, is no longer 
linear in the density as was predicted by linearization, 
cf. Eq. and cf. Eq. i|58|). but saturates toward the 
unitary limit n* m w 1. The dependence of the reflection 
coefficient R on the detuning A is plotted in Fig. I12f c1 
for Q = J and W = and in Fig. H^c) for Q = J and 
W = 2J. In the limit of a very dilute Bose-gas (see 
dashed lines for no = 0.05) we obtain the single-particle 
result T = |i(fco)| 2 given by Eq. I|17|) . showing a sym- 
metric Fano-profile for W — and an asymmetric Fano- 
profile for W — 2 J, (see also Fig. EJb,d)). We notice 
that at such weak-driving as Q = J the peak (and the 
asymmetry) in the Fano-profiles are suppressed with in- 
creasing density. However, for strong driving f2 3> 4Jno 
we recover the features of complete reflection (complete 
transmission) through the impurity site as was already 
predicted by the linearization of the impurity. 
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Hard-core Bosons 



We now consider the limit of a strongly interacting 
Bose-gas. Its Hamiltonian is given by 



H 



h.c 



U 



+A|M)(Af| + fi (\M){Q\a + a\\Q){M\ j 
+W Q al\Q)(Q\a + W M aJa |M)(M| (63) 

where the onsite-shift for two-bosons A on the same site, 
U, by far exceeds the tunneling rate J, i.e. U^S>J. Since 
double occupation of a site by two atoms A is strongly 
suppressed, we may eliminate those excitation from H, 
e.g. by imposing a? = 0. In the following we will focus 
on the limiting case U/J — ► oo, i.e. that of a Tonks gas. 
In this limit we may fermionize the Hamiltonian i|63|) 
via a Jordan- Wigner transformation (JWT) j2|j, which 
maps the commuting fields for the hard-core bosons, a.j, 
and the pseudo-spin of the impurity, er~ = |M)(Q|, onto 



anticommuting fields, 
is given by 



and /, respectively. The JWT 



JJ (l - 2c\ Cl ) , a] = cj [] (l 2^ci) , (64a) 



Kj 



Kj 



°~ = f n f 1 - 2c ^o . ct+ = / f n f 1 - 2c ^o -( 64b ) 



JV 



i*(o)> = n 



— ^ sin (k n Xj)cl 



j<0 



vac) C F 



iV 



iQ)E(- i ) s(j) n 

j<0 



n=l 



sin (k n Xj n ) a 



\vac)A, 
(66) 



where k n = nir/(M + l)a are the quasi-momenta of the 
fermionic excitations, j = (ji,j 2 , ---Jn) (with j p ^ j q ) 
denotes the position of the N bosons and (-l) s (j) the 
permutational sign of j, i.e. S(j) — J2j p >j q !■ Due to 
the cumbersomeness of the many-body wavefunction in 
terms of the bosonic operators aj , it is preferable to deal 
within the fermionic picture and extract the quantities 
of interest from the correlations for the fermions. The 
density n(xj, t) of the hardcore bosons A corresponds 
to the density for the fermions c, while the correlations, 
p(%i, Xj, t), and the momentum distribution of the Tonks 
gas, n(k,t), differ from those of a Fermi-gas, as 



n(xj, t) 

p(Xi,Xj,t) 

n{k, t) 



(atcij) = (ctcj), 

- (4 n [i - 2cj Ci 



2tt 



a x ^ 



( a m a n) 
n 



(67a) 
(67b) 
(67c) 



Q 



The fields Cj and / describe fermionic excitations for 
the new joint vacuum state of the system, |vac)c_F = 
| Q) | vac). We rewrite the Hamiltonian (|63|l in terms of 
the fermionic excitations, c, and /, and obtain 



H 



h.c. 



j<M 

+A/t/ + n(-i)^ (/t Co + h. c . 

+W Q fpclc + WmP/cIcq, 



(65) 



= X)j>o a j a j i s the number of 



where N R = Ej_>o c ] c j = Ej>o a ] 
particles to the right of the impurity site. The Hamilto- 
nian for the fermionic excitations, Cj and /, is the same 
as the one obtained for the Fermi-gas, cf. Eq. Q4U|). ex- 
cept for the appearance of the phase-factor (— 1)^" for 
the coupling Q to the impurity. 

We proceed by detailing the time-dependent scattering 
of a tonks gas with N atoms A off the impurity atom Q. 
We assume that at time t = the atoms A are trapped 
within a box of M sites to the left on impurity site. This 
corresponds to a fermi-sea of the fermionic modes Cj , and 
the state of the system at t = is given by (see also 



Diagonalizing the single-particle density matrix 
p(xi,Xj,t) one obtain the condensate fraction No(t) as 
the largest eigenvalue and the wavefunction of the quasi- 
condensate ipo(xj,t) as the corresponding eigenmode. In 
the following we denote density of the quasi-condensate 
as iio{xj,t) and its momentum distribution as rioik, t). 

In Fig. we plot the initial density for a Tonks- 
gas trapped on M — 50 sites for various filling factors 
v = 1/4, 1/2,3/4, 1, i.e. we have N = 13,25,38,50 par- 
ticles for Fig. I13f a.b.c.d). respectively. The solid lines in 
the plot on the left shows the density in position space, 
n(xj), and the dotted lines show the contribution of the 
largest eigenmode of the single-particle-density matrix 
p(Xi,Xj), no(xj). To the right we plot the correspond- 
ing quasi-momentum distributions of the gas, n(k) (solid 
line), and for the largest eigenmode of p(xi,Xj), no(fc) 
(dotted line). While the density merely resembles that 
of a homogeneous Fermi-gas with local filling factor z/, the 
momentum distributions strongly differs from the typical 
Fermi-sea (cf. Fig.[HJ) as it shows a sharp peak at k — 0, as 
one would expect from a condensate. However, the con- 
densed fraction Nq = no [xj ) (see dashed lines within 
the same figures) is not macroscopic, as it amounts only 
to sa \ZJV particles in the gas and thus the behavior of the 
Tonks-gas significantly differs from that of a true BEC. 
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FIG. 13: Initial density distribution of a gas of hard core 
bosons. The gas is trapped on M = 50 to the left of the 
impurity site, Xj = 0, and the fillings factor is (a) v = 1/4, 
(b) i/ = 1/2, (c) v — 3/4 and (d) v = 1, respectively. Shown 
are the atomic density n (solid line) and the mode of the 
quasi-condensate no (dotted line), in position (left column) 
and momentum space (right column). 



In fact we notice that the momentum distribution shows 
considerable wings, which account for the depletion in 
the quasi-condensate. Moreover, with increasing density 
v the number of particle in the quasi-condensate depletes 
considerably until vanishing completely for v = 1 (see 
Fig- 113r d)). where the system attains a Mott-insulator 
with no phase-correlations. 

We now detail the free (fi = Wq = 0) evolution of the 
Tonks gas after having opened the switch at t = 0. From 
Eq. (|65|l we obtain the state of the system at time t to 
be given by (cf. Sec II V A| l 

N 

l*t> = IQ>II [ a j( fe n,i)cj] |vac), (68a) 
n=l 

aj(k n ,t) = J — Uj tj i(t)sm(k n Xj), (68b) 

j'<0 

where Ujji (t) denotes the free single-particle propa- 
gator, cf. Eq. igBJ with W = Q = 0. The density 
n(xj,t) of the hardcore-bosons A corresponds to that of 
the fermions, given by Eq. (|43a|l . In Fig. [21 we show 
to the left the densities n(x, t) of the Tonks-gas (solid 
lines) and of the condensate mode no (xj , t) (dashed line) 
at time t = M/J = 50/ J after opening the switch. To 
the right we plot the corresponding momentum distri- 
butions, n(k,t) and no(k,t). As in Fig. ^| the sub- 



FIG. 14: Free evolution of a gas of hard core bosons. The gas 
was initially trapped on M = 50 to the left of the impurity 
site, Xj = 0, and the fillings factor is (a) v = 1/4, (b) v — 1/2, 
(c) v = 3/4 and (d) v = 1, respectively. Shown are the atomic 
density n (solid line) and the density of the quasi-condensate 
no (dotted line), in position (left column) and momentum 
space (right column) at a time t = M/J after opening the 
mirror. 



plots Fig. I14f a.b.c.d) correspond to an initial filling fac- 
tor v — 1/4,1/2,3/4,1, respectively. The density of 
the Tonks-gas, n(xj,t), corresponds to the one obtained 
for a Fermi-gas (see e.g. Fig. [Sfa) for v = 1/2), and 
shows the spreading of the gas through the impurity. 
The corresponding momentum distribution of the Tonks- 
gas, n(k,t), is shifted away from k — towards k > 
and spread in momentum space, due to the tunneling of 
the particles through the impurity site. Moreover, after 
a brief transient period starts developing an new addi- 
tional peak at k w kpjl — wit/ '2, which gains in mag- 
nitude until reaching its maximum value at t — M/J. 
This corresponds to the dynamical formation of a quasi- 
condensate uq (k, t) 2?J , which now propagates as a wave- 
packet through the impurity with k « fvr/2 (see dotted 
lines in Fig. 114(1 . The number of particle is the quasi- 
condensate, No(t) = Y^, j n o (xj , t) , gives the magnitude 
of the peak and grows with progressing time until it sat- 
urates at t w M/J. At this time the Tonks gas on the 
right side of the impurity is significantly depleted and 
one has No(t) « y/N. The dynamical formation of a 
quasi-condensate propagating with k > is clearly seen 
for filling factors v > 1/2, where it momentum peak at 
k = kp/2 exceeds the magnitude of the initial k — co- 
herent fraction. For the Mott-insulating state (commen- 
surate filling v = 1) the situation is particularly clear: 
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FIG. 15: Evolution of a gas of hard-core bosons with filling 
factor v = 1/2 for fi = and (a) W = 0, (b) W = J, (c) 
W — 2 J. Shown are the density n(xj,t) (left column), the 
momentum distribution n(k, t) and the density of the quasi- 
condensate no (xj ,t). 



purity (see Fig. I15f a)). For W — J the tunneling of 
particles through the impurity site is partially blocked 
and the dynamical formation of a monochromatic mode 
with k > is suppressed, as is see from the. Thereby 
the condensate mode no(xj, t) remains mainly localize to 
the left side of the impurity. For W — 2 J we see that 
only a small fraction of atoms passes through the im- 
purity and the momentum distribution remains centered 
at k = 0, although becoming slighty broader. We also 
see that the condensed fraction of atoms in the conden- 
sate is efficiently hindered from passing the impurity and 
remains localized to the left of the impurity site. 

In the presence of a finite coupling £1 > 0, the Hamilto- 
nian is no more bilinear in cj and Cj , due the appearance 

of the nonlinear factor (—1) Nr = Ylj >0 (l — 2ctcj), which 
makes a description of the time-dependent scattering in 
terms of the fermionic modes cj in the general case as dif- 
ficult as integrating out the full many-body Schrodinger 
equation l|63l) for the hard-core bosons A. However, the 
contribution of the nonlinear factor (—1) Nr to the dy- 
namics is negligible for strong driving, SI ^> J, |A|,W, 
and also for low densities, v <C 1. In this case the num- 
ber of atoms on the right Nn(t) -C 1, and those we set 
Nft — ► in Eq. (|65|l . The state of the system is given by 



the mott-insulator melts through the impurity (see left 
side of Fig. ITlT d)) and forms a quasi-condensate propa- 
gating as a coherent wavepacket with velocity Vk = 2 J a 
through the lattice. Thereby the initially flat develops a 
sharp peak at k = ir/a (see right side of Fig. ll4f d)). The 
number of particles in the quasi-condensate increases un- 
til reaching ~ \/N for t — Mj J. It saturates as the initial 
Mott-insulating state is significantly depleted. Moreover, 
from the evolution of the density distributions in position 
space (see plots on the right of Fig. H"4T a-d)) we see that, 
although the quasi-condensate only amounts to a fraction 
of atoms in the system, its profile no(xj,t) describes the 
transmitted part of the Tonks-gas through the impurity 
accurately, i.e. up to a multiplicative factor. 

In the presence of a finite on-site shift, Wq ^ 0, but 
with driving Q = 0, the Hamiltonian is still bilinear in 
the fermionic modes c\ and Cj . The evolution of the sys- 
tem is given by Eq. I|68[). but now with j(t) being the 
single-particle propagator for W ^ 0. In Fig. I15f a.b.c) 
we show the evolution of a Tonks gas with v = 1/2 for 
W = 0, W = J and W = 2J, respectively. On the 
left side we plot the density distribution, n(xj,t), in the 
middle the momentum distribution, n(k,t), and to the 
right the quasi-condensate no(xj,t). Darker regions cor- 
respond to regions of higher density. For no interaction, 
W — 0, the particles freely tunnel through the impurity. 
The momentum distribution is shifted towards k > as 
well as being broadened. At t « M/2J = 25/M an ad- 
ditional peak in the momentum distribution is formed 
at k w 7r/3a, which corresponds to the mode of the 
quasi-condensate, hq (xj , t) , tunneling through the im- 



N 

l*t) = H[aj(k n ,t)ct+f3(k n ,t)f^\va,c)(69 



ai(k„,t) 



■ ^2 Uij(t) sm(k n Xj), 



(69b) 



j<0 



^ UMjjt) sin(k n Xj), (69c) 



j<0 



where U a j(t) denotes the full single-particle propagator, 
cf. Eq. (|38fl . In this regime the density distribution of the 
Tonks-gas, n(xj,t), corresponds to that of a Fermi-gas 
(see e.g. the left side in Fig-JSJ. 

For the general case of many bosons A, arbitrary cou- 
pling strengths, and even finite U, we refer to the exact 
numerical simulation given in Ref. (24|. These simula- 
tions allow to test the behavior of the gas for essentially 
arbitrary repulsion U and density is, i.e. for the full 
crossover regime from a weakly interacting dilute Bosc- 
gas up to a dense Tonks gas. 



V. CONCLUSION 

We have studied a scheme utilizing quantum interfer- 
ence to control the transport of atoms in a ID optical 
lattice by a single impurity atom. The two internal state 
represent a qubit (spin-1/2), which in one spin state is 
perfectly transparent to the lattice gas, and in the other 
spin state acts as a single atom mirror, confining the lat- 
tice gas. This allows to "amplify" the state of the qubit, 
and provides a single-shot quantum non-demolition mea- 
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surement of the state of the qubit. We have derived ex- 
act analytical expression for the scattering of a single 
atom by the impurity, and gave approximate expressions 
for the dynamics a gas of many interacting bosonic of 
fermionic atoms. 

A numerical study of this dynamics based on time- 
dependent DMRG techniques, which complements the 
present discussion, will be presented in Ref. |24j . 
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APPENDIX A: SCATTERING OF GAUSSIAN 
WAVEPACKETS 

We consider the dynamics of a bosonic A-particle state 
of the form 



l*t) 



*N\ 



-i N 



vac) 



where UM,j(t) = (M\e \j) is the single particle propa- 
gator. We introduce the peak of the initial atomic density 
n = n(x ,t = 0) = NMq « N (2/tt) 1/2 5k a. 

For the scattering off the particles, i.e. for | a^o | 3* r B, 
we may neglect the finite range of the bound-state and 
have 



n m (t) 



N 



E 



a 

2^ 



= N (— 

Thus we have 

n m (t) = n Q \f t \ 2 , 
1 



Ek - A 



ft 



E k -A 



-i(k-ko)xo 



For the Fourier transform we use a saddle-point 
method, however, we have to distinguish which one is 
narrower, either the width of the wave-packet, Sko, or 
the width of dressing profile. 



2. Broad Fano- profile 



where at t — the gaussian wave-packet is given by f3o 
and 



(Xj) = J^ oe -^l{x J -x„)+ l k x^ 



where xq <C (fco > 0) is the mean position (momen- 
tum) and 5ko < it /a the width in momentum space 
at t = 0. For 5ko <C 7r/a we may take the continuum 

limit yV — > J dx/a and obtain A/b ~ (28k^a 2 /tt) , and 
Sxq w l/2<5fco 3> a/7r. The momentum representation of 
the is given by 



(k) = (a/27r) 1/2 ^, 



(2nSk 



2N-1/4 (k-k ) 2 /iSkl-i{k-k a )xa 



For a broad resonance, i.e. fltk/(Ek — A) slowly vary- 
ing on the Bloch band we expand the integral around 
k w fco, i.e. with 

1 2 

E k ~E + v (fc - fco) + 7; (k - fco) , 

2mo 



£ fc -A S -A' 

where energy £ = £fc |fc=fc , velocity v = dE k /dk \ k= k 
and effective mass mo = l/(d 2 Ek/d 2 k) \k=k - Notice 
that ink — —l/a 2 Ek for the explicit shape of the Bloch- 
band and choice of the origin in the band-middle. 
Thereby we obtain 



,,s 5k t 
dk 



-2Sk t x' t 



1. Molecular density 

For the evolution we are interested at the population 
of the molecular state, which follows as 



( (t) = (tfb) t = N 



NNq 



-Sk Q (xj —xo)+ikoXj 



where the linear propagation, spreading and the dressing 
are 



x t = x + v t, 



Skt 



5k 



1 + {25k 2 a 2 E t) 



-r, _ V 2 \t k \ 2 

VQ ~ JE—Af |fc=fc ° 



n 2 



(E - A) z + [Q 2 + W (E - A)Y a 2 /v 2 
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We notice that at i* = —xq/vq = \xq/vq\ we attain a 
maximal molecular density of 



sj 1 + (2Sk 2 a 2 E tf 



For Eq s» we might neglect the broadening/spreading. 
We recognize that T>, is maximal for such detuning A and 
initial momentum fcp where 



W 



n 2 



W 2 +v 2/ a 2 



is on the bloch-band. This corresponds to the position of 
the Fano-profile. As v 2 /a 2 = 4 J 2 — E 2 the Equation for 
E+ is implicitly cubic gives the same recursive equation. 
Thus the maximal density for a broad resonance fl ^> 
| J — A | is suppressed as 



n 2 



while far off resonance I A I 3> fi we have 



n m {U) « n Q 



(n/A) 2 



1 + (aW/v ) 2 



(Al) 



(A2) 



3. Narrow Fano-profile 

In the second case, that of a narrow resonance, i.e. a 
sharp Fano profile, we have that the dressing factor is 
narrower than the Gaussian wavepacket, hence we ap- 
proximate via a Saddle-point method. We expand the 
dressing function, 



Qtk 
E k -A 



iaW 

Vk 



-i -l 



around the momentum fc* where \T>k\ is maximal, 



r +iCi(k-k*)a- ±C 2 (k~k„) 2 a 2 



(A3) 



with 



C 
Ci 

c 2 



-in 

7* + iW im* ( .E* — A 
t I J -!- / 

7* + iW 
C 



7* V C 



7* 



c 



7* z 



c 



ra 2 + — A) (pf — jTw) 

7* 



where E 1 * = 7* = dE(k)/dka \k=k* = v{k*)/a, 

to* = d 2 E(k)/d(ka) 2 = —E(k+) are the lowest expansion 



coefficient of the dispersion relation E(k) — — 2J cos(ka) 
at k —— k+ . 

The position of its maximum fc* is obtained from the 
expansion l|A3() bv requiring 3[Ci] = with 



= 7* 



n 2 w + (E+ - A) ( 7 2 + w 2 



ft 4 + 2 (£* - A) ft 2 PF + - A) ( 7 2 + VF 2 ) 



m* 



[f! 2 + (i^-A) W\ 



7* ft 4 + 2 (£* - A) fi 2 FF + (-E 1 * - A) ( 7 2 + W 2 ) ' 

In the limit of interest (i.e. near the middle of the Bloch- 
band), we have \rrik/"fk\ = | cot(fca)| <C 1 thus we obtain 
the energy of the maximum as a series 



wn 2 



7 2 + W 2 1- + w2 y 



+ G 2 (m.). (A4) 



Since 7*, m+ = — E+ all depend on fc* the series gives 
implicitly the value of fc*. Moreover, we notice that the 
truncation to first order in to* yields the exact result for 
W = 0. 
Then with 



c 

Co 

6*2 



n 2 



7* - iW 



«7* 



7* \W 2 + -i 2 



— (w + z 7 *) /n 



1 - im^n 2 / ( 7 2 + W 2 Y 
7 2 + W 2 m+W 



7 2 + W 2 

m* 7 2 - VF 2 - 4z7*VF 



ft 2 

7 2 + W ^2\2 2 u-2 



(A5) 



o 2 y 7 * n 2 

we can compute the Fourier-integral and obtain 

l3t\A(t)(28k x(t)+i^^-) 2 + ( 



|/ t | 2 = \C \ 2 \A(t)\e 
A(t) = [l + 2(5^a 2 (C 2 +im^)] _1 , 
x(t) = xo + v+t — aCi. 



(A6) 



(A7) 



For E* w 0, i.e. near the middle of the Bloch-band, 
we might neglect the spreading of the wave-packet, i.e. 
to* fa 0. Then wc have 



I Co I 



61 



TF 2 + v 2 Ja 2 



n 2 



x 6 2 , 



C\ exp 



l+2Skla 2 C'f 



(A8) 



1 + 26k 2 a 2 C 2 
In the limit — ^ + the probability |/(£)| 2 vanishes as 



max t |/(i)| 5 



n 2 



exp 



25kl 



W 2 + v 2 Ja 2 26k 2 a 2 



(A9) 
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APPENDIX B: VARIATIONAL ANSATZ „9a,- M ( 

i\CM\ Q t = \cm\ \-J/, Uj±i,M + WSj a jt M 

For the Ansatz (|60(l we obtain the action (|61|l as 



5 = \ E - A<5 CTM |c CT | 2 + iV CT |c CT | 2 2 aj 



A - A* A 

-Oij,M + —Olj,Q, 



x V~Qj~ + Ja j+i,* + Jaj-i,<T-W8 j0 a jta ]\ A - A* N-1X 

I Y 

-nVNc* u c Q a j,M^j,Q a jt Q + c.c 



Qt -r.,« 3+liCr -r r,» }U «, lff ,j A - A* iV-lA* X*6 j0 



where the overlap s and the effective coupling A are given 
by 



with C(j = c a {t), a.j >a = aj_ (t) and Nq = Nm + 1 = N. 

Minimizing the action 5 with respect to c* and Oj, CT , 
we obtain after some algebra, s = ^^ a* M otj q, 



j 



tC * M ^W = &\cm\ 2 + X-(N-1)^±^, x = nVNc* M CQS N - l a , 

■ * dc Q x* Ar A + A* 

i c r>—^— = A — Jy — - — , 
Q dt 2 ' 
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